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Abstract 

We investigate inviscid instability in an electrically conducting fluid af- 
fected by a parallel magnetic field. The case of low magnetic Reynolds num- 
ber in Poiseuille flow is considered. When the magnetic field is sufficiently 
strong, for a flow with low hydrodynamic Reynolds number, it is already 
known that the neutral disturbances are three-dimensional. Our investiga- 
tion shows that at high hydrodynamic Reynolds number(inviscid flow), the 
effect of the strength of the magnetic field on the fastest growing perturba- 
tions is limited to a decrease of their oblique angle i.e. angle between the 
direction of the wave propagation and the basic flow. The waveform re- 
mains unchanged. The detailed analysis of the linear instability provided by 
the eigenvalue problem shows that the magnetic field has a stabilizing effect 
on the electrically conducting fluid flow. We find also that at least, the un- 
stability appears if the main flow possesses an inflexion point with a suitable 
condition between the velocity of the basic flow and the complex stability 
parameter according to Rayleigh's inflexion point theorem. 

1 Introduction 

We consider the instability in a shear flow of an incompressible viscous electri- 
cally conducting fluid with the initial velocity profile[[Q 

U = (U(z), 0, 0), U{z) ^ at z = {z u z 2 ) = (-1, +1). (1) 
*movins2008 @yahoo.fr 
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Next, we impose throughout the flow a uniform time-independent magnetic field 
B = (B, 0, 0) in the streamwise direction. The magnetic Reynolds number will 
be assumed to be small[2| i.e. 

Re m = ^y- « 1, (2) 

where L is the length scale and will be taken as the initial vorticity thickness of 
the layer, U represents the velocity scale for the flow and A = ^ stands for the 
magnetic diffusivity in which a is the electrical conductivity of the fluid and /i 
the magnetic permeability of a vacuum. 

The condition given by formula (2) is widely obtained in industrial flows or 
liquid metals, molten oxides etc... This allows one to apply the low-Z?e m ap- 
proximation (Davidson 2001) in which only the imposed magnetic field B in the 
Lorentz force expression is taken into account. This leads to the following non- 
dimensional equations 

| + (y.V)v = -Vp + ^ V 2 v + 7V(j A B), 

(3) 

V.v = 0. 

The electric current j is given byll3ll 

j = _V0 + vAB, (4) 

where <\> is the electric potential which is a solution of the Poisson equation 

V 2 = B.(VAv). (5) 

The two non-dimensional parameters appearing in eq.3 are defined as 

U L aB 2 L 

Re = , N = —— (6) 

v pU 

(Reynolds number and magnetic interaction parameter respectively). 

The magnitude of N gives information on the ratio between the Lorentz and 
inertia forces leading to the evaluation of the potential of the magnetic field which 
suppresses and transformes the perturbations. 

There are no electric or Lorentz forces generated in the non-perturbed basic 
flow @. 
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2 Governing equations 



By using the stability analysis of a shear velocity profile in the presence of a 
parallel magnetic field performed by Michael (1953), Stuart (1954) and Drazin 
(1960), we can use the normal modes for the fluctuating part of the velocity in the 
form 

v' (x, y, z, t) = \(z)exp[i(k x x + k y y - At)] (7) 

in the standard way of linear stability analysis with k x and k y representing the 
real wavenumbers in the x- and y-directions, A = u + i/3 is the complex stability 
parameter, where /3 is the growth rate of the instability and u is the frequency. 
If P > , the disturbance grows and the system becomes unstable. Whereas, if 
/3 < 0, the disturbance decays and the system becomes stable. f3 = corresponds 
to neutral instability. 

At this point, one should point out that the magnetic field stabilizes the flow 
because of the Joule dissipation action which suppresses the growing perturba- 
tions (H. But it has to be verified by an eigenvalue problem where A = u + i/3 
will be the eigenvalue. 

We consider here only two-dimensional disturbances [5] with k y = 0. For an 
arbitrary mode with k y ^ 0, the classical generalized Orr-Sommerfeld equation 
becomes 

(k x U - \){y" z - k 2 v z ) - k x v z U" + ik 2 x Nv z = -^-(v" g " - 2k 2 v" z + k 4 v z ), (8) 

Re 

for which v z = v' z = if z = (z\, z 2 ). 

As usual, k = {k 2 x + k 2 ) 2, the primes stand for the first derivatives with respect 
to z. 9 = cos~ 1 (^) is called the oblique angle between the direction of the wave 
propagation and the basic flow. 

By rearranging (8), we obtain 

(U - \){v" z - k 2 v z ) - v z U" + ikNv z = l -^(v' z " - 2k 2 v" z + k 4 v z ), (9) 

kRe 

with the following boundary conditions: v z — v' z — if z — (zx, Z2). 

We redefine new non-dimensional parameters as follow: Re = (%-)i?e, N = 
(^)ATandA = (A). 

The solution of the problem is given as a relation between A, Re, N and k in 
the form 

F(A, k, Re, N) = 0, (10) 

for any angle 9. A particular solution can be determined for a two-dimensional 
waveforms with 9 = 0. We could get a critical Reynolds number Re c which 
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corresponds to the minimum Re occurring over all k and u at which a neutral 
mode with j3 — is noticed by writing 

F c (R~e c ,N) = or Re c = G(N). (11) 

In the non-magnetic case with N = 0, the Squire theorem 11121 . |[5l requires that 
the two-dimensional perturbations are always first to become unstable since the 
smallest critical Reynolds number Re c = -jrRe c = is always for the pertur- 
bations with = 0. 

For the inviscid flow, we have 

Re -»■ oo, (12) 
and then the generalized Orr-Sommerfeld problem (9) becomes 

(k x U - \){v" z - k 2 v z ) - k x v z U" + ik 2 x Nv z = (13) 

with the condition v z = v' z — if z — (z\, z 2 ). 
By doing as (9), we have 

(U - \)(v' z - k 2 v z ) - v z U" + ikNv z = (14) 

with the condition v z = v' z = if z = (z%, z 2 ). 

Here, the non-dimensional parameters are N = (^)N, and A = (^-). 

3 Rayleigh's inflexion point for the flow 

Let us consider the linear stability of a uni-directional base flow in a channel. We 
derive the Orr-Sommerfeld equation, which governs the linear stability of uni- 
directional shear flows with respect to 3D perturbations, for viscous fluids. We 
obtain it by taking N = in (9). Then we can write 

(U - \)(v, - k 2 v z ) - v z U" = - 2k 2 v" z + k*v z ), (15) 

kRe 

with the condition v z = v' z — if z = {z\ 1 z 2 ) . 

In the inviscid case, we have the following Rayleigh's equation 

(U -~X)(v" z -k 2 v z )-v z U" =0. (16) 

Suppose that U and DU where D = are continuous in z\ < z < z 2 . Rayleigh's 
inflexion point theorem then states that a necessary (though not sufficient) con- 
dition for inviscid instability is that the base state possesses an inflexion point 
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somewhere in the domain z\ < z < z 2 . If a base state lacks an inflexion point, 
therefore, we can conclude it to be stable, for inviscids fluids. 

Consider equation (16) in the following form with the substitution v z = v, 

D 2 U 

D 2 v - (k 2 + ——=)v = 0. (17) 
U — X 

Suppose initially that the flow is unstable (/3 > 0), it is proved that an inflexion 
point i.e D 2 U = must exist for this to be so. 

Using boundary condition v(z\) = v(z 2 ) = and by making some calcula- 
tions, we get 

i 2J . r,* . d 2 u .{u ^) \ i . 1 2 , 



(k 2 + J )\v\ l dy = (18) 



where A is the complex conjugate of A. The imaginary part of this equation is 

f Z2 D 2 U\v\ 2 

From the hypothesis (3 > 0, we conclude that D 2 U must change signe somewhere 
in the domain (z±, z^). 

Then, a necessary condition for inviscid instability is the presence of an inflex- 
ion point; the absence of an inflexion point necessarily confers (inviscid) stability. 

Let us investigate what happens if the flow is affected by a parallel magnetic 
field i.e the case N ^ 0. 

The same calculation leads to 



22 

2, 



Dvrdz 



fc2 (D 2 U-ikN)(U-\) 



\U - A 



\v\ 2 dy = 0. (20) 



The imaginary part of this equation is 

" 2 U_ 

\u-\\ 2 

(fT9l) and (12TT) lead to same properties if 



/3D 2 U + kN(u - U) 2 

— \v\ ay = 0. (21) 



kN(Q - U) = (22) 

with A = Gj + i(3. 

So the inviscid flow will be unstable if there is the presence of an inflexion 
point in the main flow with the condition 

u = U(z) i.e 5R e (A) = U(z). (23) 
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4 Linear stability analysis 



We analyse linear stability of the basic flow (1) to normal mode (7). A Poiseuille 
flow with the basic profile 



U{z) = l-z 2 



(24) 



is considered. 

The eigenvalue problem (14) is solved numerically. The solution is found in 
a layer bounded at z = ±1 with U(±l) = 0. The results of calculations are 
presented in the following figures. 

For a fixed 9 = 0, we get figure 1 of (3 vs k in which a) shows the entire graph. 
b) and c) are the magnified versions of a). 

For sequential values of N, we get figures 2, 3 and 4 of (3 vs k for differents 9 
in which a) shows the entire graph, b) and c) are the magnified versions of a). 

We get also figure 5 of (3 vs k for differents N in which a) shows the entire 
graph, b) and c) are the magnified versions of a). 

For sequential values of N, we get figures 6, 7 and 8 of (3 vs k for differents 9 
in which a) shows the entire graph, b) and c) are the magnified versions of a). 
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Figure 1: (a) Growth rate (3 vs. wavenumber k with 9 — 0; (b) zoom of (a) to 
small values of k; (c) zoom of (b) to small values of k. 
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Figure 3: (a) Growth rate (3 vs. wavenumber k with N = 1; (b) zoom of (a) to 
small values of k; (c) zoom of (b) to small values of k. 




Figure 5: (a) Growth rate j3 vs. wavenumber k ; (b) zoom of (a) to small values of 
k; (c) zoom of (b) to small values of k. 
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Figure 6: (a) Growth rate j3 vs. wavenumber k with N = 0.5; (b) zoom of (a) to 
small values of k; (c) zoom of (b) to small values of k. 
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Figure 7: (a) Growth rate j3 vs. wavenumber k with N = 1; (b) zoom of (a) to 
small values of k; (c) zoom of (b) to small values of k. 
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Figure 8: (a) Growth rate j3 vs. wavenumber k with iV 
small values of k; (c) zoom of (b) to small values of k. 
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Figures 1-8 show the strong stabilizing effect of the magnetic field on the two- 
dimensional perturbations. 

It has to be stressed that the complete stabilization requires non-zero viscosity. 
It is shown in the inviscid two-dimensional analysis of Thess & Zikanov (2005) 
that the shear flow (1) cannot be completely stabilized by the magnetic field. There 
always exists a range of small k, where the flow is unstable. Such behaviour is 
in agreement with the intuitive pictures, according to which the rate of the Joule 
dissipation decreases with increasing wavelength in the direction of the magnetic 
field, and, thus, the perturbations become less and less sensitive to the action of 
the magnetic field as k — > 0. 

Typical dependence of (3 on 9 and k for the three-dimensional disturbances is 
shown in figures 2-4. The growth rate changes slowly with the wavenumber and 
the oblique angle. 



5 Conclusion 

In this paper, we revisited the inviscid instability of an electrically conducting 
fluid(modelled as a temporally evolving flow initially given by a Poiseuille flow 
velocity profile) subject to a parallel uniform magnetic field. The case of small 
magnetic Reynolds number was considered. We find an important condition be- 
tween the velocity of the basic flow and the complex stability parameter for which 
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the main flow, if it possesses an inflexion point, leads to unstability according to 
Rayleigh's inflexion point theorem. We provided detailed analysis of the linear 
instability of the problem in Poiseuille case in direct numerical simulations by 
resolving the corresponding eigenvalue problem. It shows us that the magnetic 
field has a stabilizing effect on the electrically conducting fluid; however, it re- 
mains stable for all possible values of the magnetic field since the wavenumber is 
non-zero. 
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